Intertwining operators for Sklyanin algebra and 
elliptic hypergeometric series 



A. Zabrodin * 
December 2010 



Abstract 

Intertwining operators for infinite-dimensional representations of the Sklyanin 
algebra with spins £ and —i — 1 are constructed using the technique of intertwin- 
ing vectors for elliptic L-operator. They are expressed in terms of elliptic hy- 
pergeometric series with operator argument. The intertwining operators obtained 
(VF-operators) serve as building blocks for the elliptic i?-matrix which intertwines 
tensor product of two L-operators taken in infinite-dimensional representations of 
the Sklyanin algebra with arbitrary spin. The Yang-Baxter equation for this R- 
matrix follows from simpler equations of the star-triangle type for the VF-operators. 
A natural graphic representation of the objects and equations involved in the con- 
struction is used. 
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1 Introduction 



Central to the theory of quantum integrable systems is quantum i?-matrix satisfying the 
celebrated Yang-Baxter equation. General i?-matrices with additive spectral parameter 
are parametrized via elliptic functions. The simplest elliptic i?-matrix is 

where 6a{z\T) are Jacobi ^-functions, aa are Pauli matrices, and ctq is the unit matrix. 
This i?- matrix is associated with the celebrated 8- vertex model solved by Baxter P, 
being the matrix of local Boltzmann weights at the vertex. The transfer matrix of this 
model is the generating function of conserved quantities for the integrable anisotropic 
(of XYZ type) spin-| chain. Integrable spin chains of XyZ-type and their higher spin 
generalizations can be solved by the generalized algebraic Bethe ansatz [21 13] . 

In lattice integrable models with elliptic i?- matrix fll.ip . the algebra of local observ- 
ables is the Sklyanin algebra [H |3] which is a special 2-parametric deformation of the 
universal enveloping algebra U{gl{2)). A concrete model is defined by fixing a particular 
representation of this algebra. Such representations can be realized by difference opera- 
tors. Similar to the s/(2)-case (models of XXX-type), the representations are labeled by 
a continuous parameter which is called spin, and for positive half-integer values of this 
parameter the operators representing the Sklyanin algebra generators are known to have 
a finite-dimensional invariant space. However, we allow the spin to take any complex 
value, so we are going to work in a general infinite-dimensional representation of the 
algebra of observables. 

Integrable spin chains of XXX-type with infinite-dimensional representations of sym- 
metry algebra at the sites were first studied in the seminal papers P, [7] in the context 
of high energy QCD, see also [S]. Later, lattice models with trigonometric i?-matrix 
(of XXZ-type) with non-compact quantum group symmetry were considered |9]. A 
representation-theoretical approach to models with elliptic i?-matrix and "non-compact" 
Sklyanin algebra symmetry is presently not available but there is no doubt that it should 
exist. 

In this paper we present a direct construction of the elliptic i?-matrix intertwining 
the tensor product of two arbitrary infinite-dimensional representations of the Sklyanin 
algebra. It can be realized as a difference operator in two variables, in general of infinite 
order, so we often call this object R-operator rather than i?-matrix. Another important 
object is a face type i?-matrix related to the i?-operator via a functional version of 
the vertex-face correspondence. The latter i?-matrix provides an elliptic analog of Qj- 
symbols. 

Our method closely follows the similar construction in the chiral Potts model [TUl [HI 
[T2] and the broken Zjy-symmetric model [TB| IT^ . It is based on the observation that 
the elementary L-operator is in fact a composite object built from simpler entities called 
"intertwining vectors" [TH [15] . Then the proof of the Yang-Baxter equation and other 
properties of the L-operator can be reduced to simple manipulations with the intertwining 
vectors using basic relations between them. Remarkably, all elements of this procedure 
have a nice graphic interpretation which makes them rather clear and greatly simplifies 
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the arguments. It provides simultaneously a very good illustration and an important 
heuristic tool. This graphical technique resembles both the one developed for the Chiral 
Potts and broken Zjy-symmetric models and the one known in the representation theory 
of g-deformed algebras, in particular in connection with g-deformation of 6j-symbols |16j . 

However, practical realization of these ideas in the infinite-dimensional setting is by 
no means obvious. Technically, it is rather different from what is customary in the 8- 
vertex model and its relatives. Our construction goes along the lines of our earlier work 
[T7] devoted to the Q-operator for spin chains with infinite-dimensional representations 
of the Sklyanin algebra at each site and extensively uses such really special functions as 
elliptic gamma-function and elliptic generalization of hypergeometric series. The theory 
of elliptic hypergeometric functions originated by Frenkel and Turaev in [18] is now an 
actively developing new branch of mathematics (see, e.g., [HI [201 EI] and references 
therein) . 

The elliptic -R-operator appears to be a composite object whose building blocks are 
operators which intertwine representations of the Sklyanin algebra with spins ^ and —C—l 
(the M^-operator) . They can be expressed through the elliptic hypergeometric series 
4(^3 with an operator argument. The kernels of the IV-operators are expressed through 
ratios of the elliptic gamma-function. These intertwining operators were found in our 
earlier paper [17] as a by-product of the general elliptic Q-operator construction. Here 
we re-derive this result with the help of the intertwining vectors using much more direct 
arguments. We also give a construction of vacuum vectors for the elliptic L-operator using 
the graphic technique and show how they are related to the kernel of the ly-operator. 

It should be remarked that similar results, in one or another form, can be found 
in the existing literature. In particular, the elliptic /^-operator has been found [22j in 
terms of operators which implement elementary permutations of parameters entering the 
RLL = LLR relation. A solution to the star-triangle equation built from ratios of the 
elliptic gamma- functions was recently suggested in [23]. Some closely related matters 
are discussed in the recent paper [21]. It seems to us that our approach may be of 
independent interest since it emphasizes the connection with the Sklyanin algebra and 
allows one to obtain more detailed results in a uniform way. 

The paper is organized as follows. Section 2 contains the necessary things related to 
the Sklyanin algebra, its realization by difference operators and representations. In sec- 
tion 3 we describe a space of discontinuous functions of special form, where the Sklyanin 
algebra acts, and which are identified with kernels of difference operators. Here we follow 
|17j . The technique of intertwining vectors developed in Section 4 is used in Section 
5 to construct operators which intertwine representations of the Sklyanin algebra with 
spins I and —I — 1. They appear to be the most important constituents of the elliptic 
i?-operator. In section 6 we show how the vacuum vectors for the L-operator constructed 
in [T7] emerge within the approach of the present paper. The construction of the elliptic 
i?-operator and related objects for arbitrary spin is presented in Section 7, where the 
Yang-Baxter and star-triangle relations are also discussed. Some concluding remarks are 
given in Section 8. Appendix A contains necessary information on the special functions 
involved in the main part of the paper. In Appendix B some details of the calculations 
with elliptic hypergeometric series are presented. 
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2 Representations of the Sklyanin algebra 



The aim of this section is to give the necessary prehminaries on representations of the 
Sklyanin algebra. We begin with a few formulas related to the quantum L-operator with 
elliptic dependence on the spectral parameter. 

The elliptic quantum L-operator is the matrix 



/ (2 A) So + ^4 (2 A) S3 ^2(2A)si + ^3(2A)s2 
L(A) = - I I (2.1) 

02(2A)si - 03(2A)S2 ^i(2A)so - ^4(2A)s3 



2 



with non- commutative matrix elements. Specifically, s^ are difference operators in a 
complex variable z: 

e,i2z) e,i2z) ^ ' 

introduced by Sklyanin [5J. Here 6a{z) = 9a{z\T) are Jacobi 6'-functions with the elliptic 
module r, Imr > 0, £ is a complex number (the spin), and 77 G C is a parameter which is 
assumed to belong to the fundamental parallelogram with vertices 0, 1, r, 1 + r, and to be 
incommensurate with l,r. Definitions and transformation properties of the ^-functions 
are listed in Appendix A. 

The four operators s^ obey the commutation relations of the Sklyanin algebra0: 

( — 1)"+"^/q,oSq,So = IiB-yS/jS^ — I^pS^Sp , 



(2.3) 

with the structure constants lab = (^a+i{0)9b+i{2ri). Here a,b = 0, . . . , 3 and {a,/3,'j} 
stands for any cyclic permutation of {1,2,3}. The relations of the Sklyanin algebra are 
equivalent to the condition that the L-operator satisfies the "i?LL = LLR" relation with 
the elliptic i?- matrix fll.ip . 

The parameter i in (12. 2p is called the spin of the representation. If necessary, we write 
Sa = s^^) or L*^^^(A) to indicate the dependence on i. When £ E 1^+, these operators have 
a finite-dimensional invariant subspace, namely, the space G^^ of even ^-functions of 
order Ai (see Appendix A). This is the representation space of the {2i + l)-dimensional 
irreducible representation (of series a)) of the Sklyanin algebra. For example, at £ = | 
the functions ^^{z), Osi^z) (hereafter we use the notation 9a{z) = 9a{z\^)) form a basis in 
0^, and the generators s^, with respect to this basis, are represented by 2 x 2 matrices 
(—i)^'^-^(9a+i{r]))~^aa. In this case, L(A) = -R(A — ^r/), where R is the 8-vertex model 
i?-matrix (11. ip . In general, the representation space of the Sklyanin algebra where the 
operators s^ act is called quantum space while the two-dimensional space in which the 
L-operator is the 2x2 matrix is called auxiliary space. 

As is proved in [25j, the space 64^ for i E |Z_|_ is annihilated by the operator 



W. - cT(-l)' \^^ + ^] ^i(2^ + 2(2f-2fc + l)r^) (2,.2.+i).a. .3 4) 



k=0 



^The standard generators of the Sklyanin algebra |4|, Sa, are related to ours as follows: Sa 
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where c is a normalization constant to be fixed below. Hereafter, we use the "elliptic 
factorial" and "elliptic binomial" notation: 

The defining property of the operator established in is that intertwines 
representations of spin £ and of spin + 1): 

W, = si'^'^^We , a = 0, . . . , 3 . (2.6) 

The same intertwining relation can be written for the quantum L-operator fl2.1l) : 

(A) = L(-^-i) (A) W^ (2.7) 

Note that the operator serves as an elliptic analog of {d/dzY^^^ in the following 
sense. In the case of the algebra s/(2), the intertwining operator between representations 
of spins £ and — £ — 1 (realized by differential operators in z) is just {d/dzY^^^. It 
annihilates the linear space of polynomials of degree < 2£ (which results in the rational 
degeneration of the elliptic space O^^). 

For us it is very important to note that can be extended to arbitrary complex 
values of C. in which case it is represented by a half-infinite series in the shift operator 
^2rid^ [17j. The series is an elliptic analog of the very- well-poised basic hypergeometric 
series with an operator argument. The explicit form is given below in this paper. The 
intertwining relations (12. 6 p hold true in this more general case, too. 

Very little is known about infinite-dimensional representations of the Sklyanin algebra. 
The difference operators (12. 2 p do provide such a representation but any characterization 
of the space of functions where they are going to act is not available at the moment, 
at least for continuous functions. On the other hand, the difference character of the 
operators (12. 2 p suggests to consider their action on a space of discontinuous functions of 
a special form. The latter are naturally identified with kernels of difference operators. 
This formalism was used in our earlier paper [T7]. It is reviewed in the next section. 



n 
m 



[n]\ 



[m] \[n — m] 



(2.5) 



3 Kernels of difference operators 

Let 6{z) be the function equal to zero everywhere but at z = 0, where it equals 1: 6{z) = 0, 
2; 7^ 0, 6{0) = 1. (We hope that the same notation as for the conventional delta-function 
will cause no confusion because the latter will not appear in what follows.) Clearly, 
z5{z) = and S\z) = 5{z). 

Consider the space C of functions of the form 

f{z) = Y.fk6{z-u + 2kr]), /, gC, (3.1) 

fcGZ 

where G C. This space is isomorphic to the direct product of C and the linear space 
of sequences {fk}kez- We call functions of the form (13. ip combs. Clearly, the Sklyanin 
algebra realized as in (12. 2p acts in this space (shifting z/ — )■ ± r/). 
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A comb is said to be finite from tfie riglit (respectively, from tfie left) if there exists 
M G Z such that = as > M (respectively, k < M). Let (respectively, C^) be 
the space of combs finite from the left (respectively, from the right). 

We define the pairing 

{F{z),6{z-a)) = F{a) (3.2) 
for any function F{z), not necessarily of the form (13. ip . In particular, 

{6{z-a), 6{z-b)) = 6{a~b). (3.3) 

Formally, this pairing can be written as an integral: 

{F{z), 6{z -a))= f dzF{z)6{z - a) (3.4) 



(perhaps a g- integral symbol would be more appropriate). We stress that the integral 
here means nothing more than another notation for the pairing, especially convenient in 
case of many variables. By linearity, the pairing can be extended to the whole space of 
combs. We note that the pairing between the spaces and is well defined since the 
sum is always finite. 

Combs are to be thought of as kernels of difference operators. By a difference operator 
in one variable we mean any expression of the form 

D = ^ Cfc(2)e(^+2'='')^^ , /i G C . (3.5) 
feez 

The comb 

D{z, = E ^kizMz - C + + 2kr]) , (3.6) 

fcez 

regarded as a function of any one of the variables z, (, is the kernel of this difference 
operator in the following sense. Using the pairing introduced above, we can write: 

(D/)(^) = / D{z, OfiOdC = E Ck{z)f{z + /X + 2kv) . (3.7) 

feez 

The kernel D{z, () can be viewed as an infinite matrix with continuously numbered rows 
(z) and columns {(). Then the convolution with respect to the second argument of the 
kernel, as in (13. 7p . defines action of the operator from the left. The convolution with 
respect to the first argument defines the action from the right, 

(/D)(^)= J f{OD{C,z)dC, (3.8) 

equivalent to the action of the transposed difference operator from the left: 

= E e-('^+2^'')^^Cfc(z) = J2 Ckiz - /i - 2A:77)e-(^+2W . (3.9) 
fcez fcgz 

The transposition t is the anti-automorphism of the algebra of difference operators such 
that (c(z)e°'^^) = e~°'^^c{z). In terms of the above pairing we can write (/, Dgf) = 

(D7,^). 
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The following simple remarks will be useful in what follows. Let F{z),G{z) be any 
functions, then F{z)D{z,()G{(), with D{z,() as above, is the kernel of the difference 
operator 

FBG = J2 Ckiz)F{z)G{z + /i + 2kri)e^^^^^'^^^^ 
fcez 

which is the composition of the multiplication by G, action of the operator D and sub- 
sequent multiplication by F. Let D^^\zX)-, D'^'^\zX) be kernels of difference operators 
13(1)^ 13(2) respectively, then the convolution 



JdCD'~'\z,OD^'\^,0 



is the kernel of the difference operator D^^^D^^^ If the kernels D^^\z,C), D^'^\z,C) are 
combs finite from the left (right) as functions of z, then the convolution is always well 
defined and the resulting kernel belongs to the same space of combs. 

The kernels of Sklyanin's operators 02.21) are: 

-<^' .) . ^..ig^ - . . - ?^it^ _ . . , . ,3.10, 

Note that Sa{—z, —z') = Sa{z,z'). Let us find the kernel of the L-operator (12. ip . Using 
identities for theta-functions, it is easy to see that 

mX) = e,{2X + 2iv)V~\X + ir^,z)(^^'-^ + '^^ 5(.-C-r]) ) ~ ^^'^^^ 



where V^(A, z) is the matrix 

V{X,z) 
and V~^{X,z) is its inverse: 



9i{z + X) 9s{z + X) 
Uz-X) U^-X) 



^ ' ' 2e^{2z)\-eA{z-x) h{z + x) ) 

(recall that da{z) = 9a{z\^)). A crucial point is that the diagonal matrix with delta- 
functions factorizes into the product of column and row vectors: 



and thus so does L^{X). The vectors which represent the f actor ized kernel of the L- 
operator are "intertwining vectors" introduced in the next section. 



4 Intertwining vectors 

We introduce the 2-component (co)vector 

10 = ( If^l ) , (CI = (^"4(0, ^^3(0). (4.1) 
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Figure 1: Intertwining vectors. 

The vector orthogonal to (0 is IC)"*" = ( ^^/^l h covector orthogonal to \0 is 
HC\ = -e.iO), so (CIC)^ = ^(CIC) = 0. More generally, we have: 

(eic)^ = 2^i(e+c)^i(e-c) = -^(eic)- 

Note also that 



|C+i(l + r)) = e-^-2-^|C)^ 
Introduce now the intertwining vectors 



2^1 (2;^) 



+ X)S{z - z' + rj) + \z - X) S{z -z! - rj)), 



, ^ i\z - X)^5{z -z' + r])-\z + X)^5{z - z! - r))) 
^J2e,{2z)^ ' 



and the corresponding covectors 



>:'(A) 

It is easy to check that 



-jl—Uz + X\b{z -z'-rri)-r{z-X\ b{z - z' - 7?)), 

Y 2^1(22;) 

(Hz - X\5iz -Z' + Tj)- Hz + X\S(Z -z! - Tj)). 

2^1(2^)^ ^ 



ei{2z') 
\ 0i{2z) 



eii2z') 



K (A + V)) , 



(4.2) 
(4.3) 



(4.4) 
(4.5) 

(4.6) 
(4.7) 
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Figure 2: The graphic representation of the relation W^'^{\ — iJ,)(^(f)l,(X + |) — 

= W^'''''{X — ijL)(^(f)l,(iJi + I) 0^/(A — |)^. The horizontal bold line segment common 
for the covector to the left and the vector to the right means taking scalar product of 
the two-dimensional (co)vectors. The intersection point of the spectral parameter lines 
corresponds to the "vertex" W^''^{\ — /i). 



The intertwining vectors satisfy the following orthogonahty relations: 

[<Pl,{X)\rA>^)) = e^{2X)6{z' - z"){6{z-z'+v) + S{z-z'-7i)), 

ei{2z) 
'0i{2z') 

' 1 
1 



'(A + 77)1 C(A - y))) = ^i(2A)^^ S{z' - z") {S{z-z'+r)) + S{z-z'-r) 

ei{2X) 



dC 



dC 0c(A))(0^(A) 
0i{2C) 



0i{2z) 



1 
1 



(4.8) 
(4.9) 

(4.10) 

(4.11) 



The general scalar product of two intertwining vectors is 



0^'(A) ^Ui^) 



^4^i(2z)^i(2C) 

X [0,{z+C+X-fi)0,{z-C+X+fi)5iz-z'+r))S{C-C+v) 
-e,{z+C+>^+l^)0i{z-C+X-i^)5{z-z'+v)6{C-C-v) 
+ e^{z+C-\-lJ^)0l{z-C-X+^i)6{z-z'-r])6{C-C+l^) 

-e,{z+C-X+pi)9,{z-C-X-pi)5{z-z'-rj)S{C-C-v)}- 

It is a matter of direct verification to see that such scalar products satisfy the "intertwin- 
ing relation" : 

W^'^iX - n){<Pl,{X + |)| 4(/. - |)) = W^''^'{X - n){<Pl,{^i + |)| 4(A - |)) , (4.12) 
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1^ 



Figure 3: The kernel of the L-operator L5(A,/i) = (pfi^ 



where the quantities ly^'^(A) solve the following difference equations in z,(- 



ei{z + ( + \ + r]) 



^ ^ ei{z + c + x + 7]) ^ ^ 



(4.13) 



These equations can be solved in terms of the elliptic gamma- function r(z|r, 2ri) := {'{z) 
[261 127] (see Appendix A): 



W'^iX) 



-2.ix.ln r(2: + C + A + 77)r(2:-C + A + ^7) 
r(2 + C- A + r/)r(2;-C-A + r/) ' 



(4.14) 



There is a freedom to multiply the solution by an arbitrary 2?7-periodic function of 2 + C 
and 2 — C- We put this function equal to 1. (However, this does not mean that this is the 
best normalization; other possibilities will be discussed elsewhere.) In our normalization 



iy"'^(A)iy"'^(-A) = 1 



(4.15) 



but ly^'^(A) is not symmetric under permutation of z and C^. 

The intertwining vectors can be represented graphically as shown in Fig. [H The 
vertical line carries the spectral parameter and serves as a line of demarcation between 
the "real" (transparent) world and the "shadow" world. Then the relation fl4.12p means 
that the horizontal line in Fig. [2] can be moved through the intersection point of the 
two spectral parameter lines. This intersection point is a new graphic element which 
corresponds to W^'^^iX — /i). 

The kernel of the L-operator for the representation of spin i can be written in the 
factorized form as the product of intertwining vectors: 



A± = A±(£ + ^)r7. (4.16) 



It clear that the spectral parameter A and the representation parameter Iri enter here 
on equal footing, so the notation L'^^)^(A) = L^(A+,A_) is sometimes also convenient. 
Graphically, the kernel of the L-operator is shown in Fig. |3l 
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Figure 4: The vertices W^'''{X — fi) and W^{X — fi). 

Intertwining operators for arbitrary spin 



There is a relation which is "dual" to ( I4.12p (see also Fig. ^ meaning that it can be read 
from the same configuration of lines in the figure by exchanging the real and shadow world 
pieces of the plane (see Fig. [5]). Two new elements appear: first, the vertex W^{X — fi) 
is different from the one in Fig. [2] and, second, one should take convolution (/ d() with 
respect to the "intermediate" variable ( associated to the finite triangle in the shadow 
world. The two vertices, W^'''{X — fi) and W^{X — fi), are shown separately in Fig. HJ 
According to Fig. [5l the dual relation has the form 

/rfCl^C^(A-/i)|0S(A-|))(0S(/^+|)|=/rfCW^5(A-/i)|0^(/x-|))(0^(A + |)|. (5.1) 

Changing the notation A — t- A+, /i — t- A_, one can write it as 

J dC W^%X+ - A_)4(A+, A_) =JdC W^^X^ - A_)L^(A_, A+) 

which is just the intertwining relation for the L-operator L''^^(A) = L(A+, A_) (12. 7p . with 
W^{X^ — A_) being the kernel of the difference operator W^. Taking this into account, 
we are going to find solutions for the in the space of combs finite either from the 
right or from the left. 

Let us take the scalar product of both sides of equation (15. ip with the covector 



0^(A + f^ 



from the left and the vector 



onality relations (14. 8p . (14. 9p . we obtain: 



4>l'iX + ^)) from the right. Using the orthog- 



This functional relation for can be solved in terms of W^''^ with the help of (I4.12p : 
W^{X)W'''^{X + r]) = 9i{2(). However, this solution is not exactly what we need because 
it is not a comb-like function. Proceeding in a slightly different way, one can rewrite 
(O system of difference equations for W^{X): 



w^+^(x\ - ^i(2C + 2r/) ^i(^ + C-A) 



^i(2C) ^i(z + C + A + 2r/) 
e^{2C + 2r,) ^i(C-^-A) 



W^c(A), 



(5.3) 



wnx). 
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Figure 5: The graphical representation of equation Ii5.1\) . 



Comparing with f l4.13p . one immediately finds a solution in the space of combs: 

with VT^'^ given by f l4.14p and arbitrary z/. (The factor in front of the sum is also a solution 
but in the space of meromorphic functions.) The function c(A) introduced here for the 
proper normalization is not determined from the difference equations. It will be fixed 
below. One may truncate the comb from the left choosing z/ = A; then the coefficients in 
front of 6{z — ( — X + 2kri) with A; < vanish because the function W'''^{X + 77) has poles 
a.t ( = z — X + 2kr], k < —1. Another possibility is to truncate the comb from the right 
choosing u = —A; then the arguments of the delta-functions at k < exactly coincide 
with the half-infinite lattice of zeros of the function W^'^{X + rj), and so one can make 
the truncated comb by taking residues. Below we use the first possibility and consider 
the solution 

(5.4) 

^ ^i(2z-2A + 4fcr7) x ^ a 

= c(A) > — — r— T — d{z - C - X + 2k'r]) 

which is the kernel of the difference operator 

= c(A) J2 e^^'^^+'i')^'-^+^'''^y'Wi {2z -2X + 4kr]) 

k>0 

^ \-{2z-2X + 2kr^)r{-2X + 2k7^) 
r{2z + 2ri + 2kri)r{2ri + 2kri) 
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Rewriting the coefficients in terms of the elhptic Pochhammer symbols with the help of 
flAlSp . f lA17p and extracting a common multiplier, we obtain 



z— A I or,! \ z—Xl r A] 



g2.iA(.-A)/„r(_2^)r(22 - 2A + 2ri) ^ + 2A;] 



w(A) = c(\)- ^ , y 



r(2r^) r(2^ + 2r^) ,to [§ + 1] 



k 



where c(A) = ie'^'o ri£){T)c{\). The infinite sum can be written in terms of the elliptic 
hypergeometric series (see Appendix A for the definition) with operator argument: 

-,,,e2-^(^-^)/''r(-2A)r(2z-2A + 2r/) fz-X X 2nd\ -xa 
W(A) = c(A) \ ^/ ^ ^ •4^3 : : e^'''^" re ^''^ (5.5) 

Here the double dots mean normal ordering such that the shift operator e"^^^^^ is moved 
to the right. By construction, this operator satisfies the intertwining relation 

W(A - /i)L(A, /i) = L(/i, A)W(A - /i). (5.6) 

The intertwining property (15. 6p suggests that W(A)W(— A) = id or, equivalently, 

JdCW^{X)W^,i-X) = 6{z-z') (5.7) 

which is a shadow world analog of (14.151) . This is indeed true provided that the function 
c(A) is fixed to be 

ax, - (5.8) 

where the constant po is 

r(2,) efi(^-'-3-) 

Po — — —T — - — jTT-r l^-yj 
ie—r]D{T) ir]D[2T]) 

(clearly, there is still a freedom to multiply c(A) by a function ^{X) such that ip{X)(p{—X) = 
1). It should be noted that the very fact that the product W(A)W(— A) is proportional to 
the identity operator is by no means obvious from the infinite series representation ( 15. 5p . 
This fact was explicitly proved in |22] with the help of the Frenkel-Turaev summation 
formula. For completeness, we present some details of this calculation in Appendix B. It 
is this calculation that allows one to find c(A) explicitly. 

We thus conclude that the properly normalized intertwining operator W(A) reads 



TTiX I 2TviXz 



r(2z-2A + 2r/) z-X X ^vd \ -xo 

or, in terms of the parameter (i = 2£ + 1 G C related to the spin C. of the representation, 

W, = W{dr]) = e--^^^+2-'^-Ei^i;:lil^iM .^,^3 (^_d; -d; e^"^^^ • e'^'^^^ (5.11) 

\{2z + 2ri) ' \r] J ' 

It is not difficult to see that the change of sign z — )■ —z transforms W(A) to another 
intertwining operator for the Sklyanin algebra, which is an infinite series in shifts in the 
opposite direction. (It is this latter operator which was constructed in the paper |17j.) 
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It can be obtained within the same approach if one uses the other possibihty to truncate 
the comb which has been discussed above. If £ G (i.e., d G Z+), then the elliptic 
hypergeometric series is terminating and both operators are represented by finite sums 
(containing d+ 1 terms). Moreover, they coincide with each other and are explicitly given 
by the formula 

W,= (.e-(-''+i)r/,,(r))'X:(-l)' 

which coincides with equation f l2.4p . 

Let us conclude this section by summarizing the graphic elements of the diagrams and 
rules of their composing. The plane is divided into "transparent" and "shadow" pieces 
by a number of straight dashed lines in such a way that each segment of any line is a 
border between pieces of the different kind. Each dashed line carries a spectral parameter 
denoted by A, /x, etc. There may be also bold straight lines which become dotted when 
they go through shadow pieces of the plane. Each shadow piece (bounded by dashed or 
dotted lines or by infinity) carries a complex variable denoted by (, etc. Those which 
sit on infinite pieces are fixed while those which sit on finite pieces bounded by lines 
of any type should be "integrated" in the sense of the pairing (13.41) . The intersection 
points of the dashed lines are of two types depending on the way how the transparent and 
shadow parts are adjacent to it. Correspondingly, there are two types of vertex functions 
shown in Fig. |H The intersection of a dashed line with a bold one corresponds to an 
intertwining (co)vector as shown in Fig. [T] Finite bold segments mean taking scalar 
products of (co)vectors associated with their endpoints. 



d 
k 



ei{2z - 2{d - 2k)r]) 
nto^i(2^ + 2(A;-j» 



(5.12) 



6 Vacuum vectors 

In order to make a closer contact with our earlier work [17], it is useful to demonstrate 
how the vacuum vectors for the L-operator can be constructed within the approach 
developed in the previous sections. Let us recall the general definition of the vacuum 
vectors. Consider an arbitrary L-operator L with two-dimensional auxiliary space C , 
i.e., an arbitrary 2x2 operator- valued matrix 




The operators Ljj act in a linear space "H which is called the quantum space of the L- 
operator. For the moment, let 0, ip, etc denote vectors from and X, Xi, etc vectors 
from Ti, then acting by the quantum L-operator on the tensor product X ® 0, we, 
generally speaking, obtain a mixed state in the quantum space: LX®0 = Xi 001-1- X2®02- 
The special case of a pure state, 

LX®0 = X'®V, (6.1) 

is of prime importance. The relation (16.11) (in the particular case "H = C^) was the 
key point for Baxter in his solution of the 8- vertex model [IJ. (This is what he called 
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the "pair-propagation through a vertex" property.) Taking the scalar product with the 
vector tp-^ orthogonal to ip, we get: 

(V^^L0)X = O, (6.2) 

i.e., the operator K = L0) (acting in the quantum space only) has a zero mode 
X eT-L. Suppose (16.11) (or (16.21) ) holds with some vectors 0, ip; then the vector X is called 
a vacuum vector of the L-operator. An algebro-geometric approach to the equation (16.11) 
for finite- dimensional matrices Lj^ was suggested by Krichever [28] and further developed 
in [29| [30] . In our paper [T7] the Baxter's method of vacuum vectors was adopted to the 
infinite-dimensional representations of the Sklyanin algebra. 

For L-operators with elliptic spectral parameter it is convenient to pass to the elliptic 
parametrization of the components of the vectors (p, ip as is given by (14. ip . Writing 
L(A)|C) (respectively, (C|L(A)) we mean that the 2x2 matrix L acts on the 2-component 
vector from the left (respectively, on the 2-component covector from the right). Similarly, 
we introduce right and left vacuum vectors Xr, Xl according to the relations 

(CIL(A)X^ = {^\X'^ , X^(C|L(A) = . (6.3) 

In the latter formula the matrix elements of L act on Xl from the right. Introducing the 
operator 

K = K(C,0 = (C|L(A)|er, (6.4) 

we can rewrite (16. 3p as KX/j = X^K. = 0. The explicit form of the operator K can be 
found from ([21]), (Q: 

K = K(C, = p{z)e^'' + p(-^)e-''^^ , (6.5) 

where 

P(^) = n + <-K + l)Oi{z + ee + A_ + I 

These difference operators appeared in [301 HI] later were independently introduced 
in [311 ES] . So, the equations for the right and left vacuum vectors read 

p{z)Xr{z + r]) = -p{-z)Xniz - r]) , (6.6) 

p{-z - r])XL{z + r]) = -p{z - 7])Xl{z - r]) . (6.7) 

Instead of solving these equations explicitly, below we show how the vacuum vectors 
emerge within the approach of the present paper. The key relation is (see Fig. |6|) 



/ dC (0|'(/i + |)|0^(A+ - |))(0^(A_ + |)| W^'^iX^ - p)W^,{X^ - /i) 
dC (0|'(A+ + ^)|4(A_ - f ))(4(/i + h\w^''%X+ - P)W^'{X- - P). 



The left-hand side represents the action of the L-operator L(^^(A) = L(A+, A_) to the 
covector ^</>|'(A+ + |) in the auxiliary space from the right and to the vector W^'^{\j^ — 
jj)W^,{X^ — /i) in the quantum space from the left. It is convenient to denote 

X^/{z\X^, A_) = W^'\X^ - |) W^^X^ - |), (6.9) 
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X X 

+ 



+ 



1^ 



Figure 6: The graphic representation of equation ^6.8\) : action of the L-operator to right 
vacuum vectors. 

then relation (16. 8p can be rewritten (after setting z' = ^ ±ri and some transformations) 
as the system of equations 



where X^/ = X|j'«'(z|A+, A_) and 





. + A_+r7)0i(e-r-A4 


_-A_ 


_-2^) 














-A_ 


-2/i) 














- + A_ 


- + 2/i) 












-A_-r/)^^i(e-f+A+ 


+ A_ 


+ 2/i) 





(6.10) 



We note that setting fi = one obtains from (I6.10p 
e|L(A+, A_)Xi'«' = (e'l (aoXr'''«'+'' + KXi^^'^'^^) = (coXr^'^'+^ + d^Xlf^^^'-^ 

where = a(/i = 0), etc. This means that X^'^ (z|A+ — |, A_ — |) is the right vacuum 
vector for the L-operator (see the first equation in (16. 3p . Moreover, we conclude that 



(6.11) 
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Figure 7: The intertwining relation RL®L = L®LR. 
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Figure 8: The kernel of the R-operator R^^',{X+, X_\^+, /i 



One can that the vacuum vector is in fact a composite object. It is a product of two 
ly-functions. Equations f l4.14p . f l5.4p together with the 3-term identity for the Jacobi 
theta-function imply the relation 

e L(A+, X^)X^/iz\X^, A_) = 0i(2A_ + X|'«'(^|A++r/, X^+r]) (6.12) 



which is equation (4.22) from our paper [T^ written in the slightly different notation. 
The left vacuum vectors can be considered in a similar way. 



7 The i?-operator and related objects 

The i?-operator R = R(A+, A_|/i+,/i_) intertwines the product of two L-operators: 
R(A+, A_|/i+, /i_)L(A+, A_) ® L(/i+, /i_) = L(/i+, /i_) ® L(A+, A_)R(A+, A_|/i+, /i_). [7.1] 
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Passing to the different notation, A± = A ± (£ + |)?7, /i± = ± (£' + ^)ri, we can rewrite 
fl7.ip in a more conventional form: 



R(^^')(A,^) lW(A) ® L(^')(/.) = L(^')(/i) ® lW(A) R(^^')(A,/i). 



(7.2) 



Here R(«')(A,/i) = R(A+,A_| is a difference operator in two variables acting in the 

tensor product of the quantum spaces for the two L-operators. In terms of the kernels 
equation (17. ip reads: 



JdC JdCRll',iX+,\-\f^+,f^-)L'liKA-)L'^^,ifx+,fi-) 
JdcJ dCmfi+, /i-)L^:(A+, A_)i?g:(A+, A_| 



(7.3) 



Graphically it is shown in Fig. [71 The figure clarifies the structure of the kernel of the 
i?-operator which is shown in more detail in Fig. El It is clear that the kernel is the 
product of four W^-vertices: two of them are of the iy^'''-type (meromorphic functions) 
and the other two are of the W^-tjpe (comb-like functions). Specifically, we can write: 

i?^^/(A+, A_|/i+, /i_) 
= W^'^'iX^ - fi^)W.^:i\- - /i-) W(A+ - /i+)W^«'(A_ - ^+) 



c(A_ -/i_)^i(2C') V- / A/ X o,/ N 



X 



cA+-/i+gi2C 
H^f'^(A+-/i+ + r/) ^ 



which is the kernel of the difference operator 

R(A+,A_|/i+,;U_) = Vr^'^'(A+ -/i_)W("')(A_-/i_)W(")(A+-^+)l^"'"'(A_ -/i+) (7.4) 

(here the notation W*^^) means that the operator W acts to the variable z. In full, the 
i?-operator reads 



(A+-/.+ )2-^(A_-M_)2 + i2li{A+-/.+)^+^{A_-/._)2' 



R(A+, A_|/x+,/i_) = e 



^ _2ii(A^_^_)^ V{z + z' + \+ - ^i- + ri)V{z - z' + X+ - fi- +r]) 
^ " riz + z' - X+ + fi_ + r])r{z - z' - X+ + +r]) 



X 



r(2/-2(A_-/i_) + 2r/) ^ (z' + - A_ _ - A_ _ 



r(2^' + 2r7) 



.41^3 



-(A_-/._)9,, 



r(2z-2(A+-/i+) + 277) . / z + /x+-A+_ fi+-X+_ 



r{2z + 27]) 



V 



V 



^ ^-izi(x—i,+)z r(^ + z' + X_- + r])r{z - z' + X_ - fi+ + T]) 
r{z + z' - X- + fi+ + r])r{z - z' - X- + fi+ + r])' 



(7.5) 
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Figure 9: The Yang-Baxter equation for the R-operator. 



The difference operators in the third and the fourth hnes of the r.h.s. commute because 
they act in different variables but both of them do not commute with the operator of 
multiphcation by the function W^'^ (A_ — Note that the /^-operator can be also 

written in terms of the ^cor^ series due to the identity 



1 4(^3 



z-\ X 



\ 1] 1] 21] 17] J 

(7.6) 

The Yang-Baxter equation for the i?-operator is schematically shown in the self-exp- 
lanatory Fig. [9l One can see that as soon as the i?-operator is a composite object, the 
Yang-Baxter equation can be reduced to simpler equations for its elementary constituents. 
The latter are the W^-vertices of the two types. For them one can prove a sort of the 
star-triangle relations 

W''^\^l - u)W''''"{\ - f^)W^„i\ -^) = I dCW^iX - n)W''^^i\ - u)W^„{n - v) (7.7) 
W'^'\\ - n)W'"'''{fi - v)Wln[\ -v)= \ dQA/n^ - v)W^'''[\ - z/)W^f„(A - /x) (7.8) 



schematically shown in Fig. [TOl The proof is given in Appendix B. As is seen from Fig. 
[HI the proof of the Yang-Baxter equation is reduced to sequential transferring of vertical 
lines from the left to the right through intersection points of the other lines with the 
use of the star-triangle relations (17. 7p and (17.81) at each step. Let us note that the both 
sides of the star-triangle relations (17. 7p and (17. 8p represent the kernels of the difference 
operators explicitly written in Appendix B (( IBSP and (]B6P respectively). 

There is an object "dual" to the i?-operator R in the sense that its kernel is graphically 
represented by the same pattern, with shadow parts of the plane being complimentary to 
those in Fig. [HI This duality provides a transformation which is an infinite-dimensional 
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X Z" V 



II Z" V 



Figure 10: The star-triangle equation \1. 7\ ) for the W -operators. Equation ( [7.g| j corre- 
sponds to the same configuration of lines with complimentary shadow parts of the plane. 

version of the vertex-face correspondence. It sends the /^-operator to a difference operator 
in one variable rather than two. We call it the S'-operator. It acts in the variable z 
and depends on z' and z" as parameters. Its kernel, S^{z', z"\X^, A_; or simply 

S^{z',z") in short, is shown in Fig. [TTl This kernel is to be regarded as an i?-matrix 
for a face-type model with complex variables associated to shadow parts of the plane. It 
generalizes the fused Boltzmann weights of the SOS-type 8- vertex model ^3\. According 
to Fig. [11] it reads 

Sliz', z") = J dCW^^iX^ - ^^-)W'''HK - f^+)W^''"{X- - f^-)W^{X- - ^+). (7.9) 

The convolution is taken with respect to the variable sitting in the finite parallelogram 
at the center of Fig. [TTl Since each of the two PF^^-vertices is represented by a half- 
infinite sum of the type (15.41) . the whole expression (17. 9p is a double sum. Performing 
the convolution and re-arranging the double sum, we can write 

S|{z',z") = c{X+-fl^)c{X^-^l+)0^{2OY.M^,^',z")6{z-^-X+-X.+f^++^_ + 2n7]), 

n>0 

where 

An{z, z' , z") 

" gi(2^-2A+ + 2/i_+4fcr7) W''''-^++^'-+^'''>{X+-^i+) iy^-^++^-+^^'''""(A- -/x.) 
Using the explicit form of the VT-functions it is straightforward to show that the kernel 
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Figure 11: The kernel S^{z',z") "dual" to the kernel of the R-operator (of. Fig. 
S^{z',z") is expressed in terms of the elliptic hypergeometric series loUg as follows: 



Sl{z',z") = C'6'i(2^)e^''*^"^^^^^+~^"^^^^^""''+^^"^^''+"^+^^'^ 

6i{z - z' + - fi+ + T]) 9i{z - ^ + yu+ + - A+ - A_) 
ei{z-z'+^i_+^i+-2\++7]) ^i(2-^-/i++/i_-A+ + A_) 

V{2z + 2/i_ - 2A+ + 2-q) V{2aji]) 



X 



(7.10) 



r(2^ + 277) n r(2(ai-a, + l)r7) 

X ioC<;9(ai; 04, ... , aio) 5{z-^- X+- \^+ fx++ n_ + 2nr]). 

n>0 

Here C is a constant which depends on the spectral parameters, 6i{x) = 6i{x\2r]) and 
the values of the aj's are 

2; + /i_ — A+ /i_ — A+ Z ± ^ + yU+ + yU_ — A+ — A_ 

ai = , 04 = , 05,6 ^ 



rj rj 2ri 

z ± z' + n_ — n+ + 7] z ± z" + X- — X+ + rj 

0^7.8 = 7: , ttg.io = 7: 

2ri 2r] 

One can see that the series iqUq is balanced (the balancing condition (]A21|) is satisfied) 



and terminating (ag = because of the ^-function). Equation f l7.10p is a version of the 
Frenkel-Turaev result [H] adopted to continuous values of parameters and obtained by 
a different method. The S-operator satisfies a sort of the Yang-Baxter equation which 
can be graphically represented like in Fig. [9] with transparent pieces of the plane being 
changed to the shadow ones and vice versa. 

Another object closely related to the i?-operator is the "transfer matrix on 1 site" 

T(A+, A_|;U+,;U_) = tr^ (^R(A+, A_|;U+,/i_)P) , (7.11) 
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where P is the permutation operator of the two quantum spaces and the trace is taken in 
the space associated with the spectral parameters fi±. The kernel of this transfer matrix 
is 

r|(A+,A_|/z+,/i_) = J rfe^g(A+,A_|/i+,/i_). (7.12) 

It is not difficult to see that this kernel is expressed through the kernel S^{z',z") given 
by (\7J0\\ as follows: 

r/(A+,A_|/i+,/i_) = S'|(e,^)(A„,A+|/i+,/i_). (7.13) 

(Note the exchange of the spectral parameters A+ O A_ in the right-hand side.) The 
easiest way to see this is to draw the corresponding pictures. 



8 Concluding remarks 

In this paper we have presented a unified approach to intertwining operators for quantum 
integrable models with elliptic i?-matrix associated with the Sklyanin algebra. We work 
in the most general setting of infinite-dimensional representations (with a complex spin 
parameter i) realized by difference operators in the space of functions of a complex 
variable z. The elementary building blocks are so-called intertwining vectors and W- 
functions which are defined in terms of their scalar products. These elements have a 
nice graphic representation as diagrams in the transparent/shadow plane which allows 
one to easily construct more complicated objects like L-operators, their vacuum vectors 
and different kinds of i?-matrices and to prove relations between them. An important 
constituent of the construction is the intertwining operator for representations with spins 
i and For general values of i, it is given by the elliptic hypergeometric series 

with operator argument. 

In fact the material presented here is only the very beginning of the theory of inte- 
grable "spin chains" with elliptic i?-matrices and infinite-dimensional space of states at 
each site. Indeed, our discussion has been focused on a single L or R operator which is 
relevant to a spin chain of just one site. The next step is to construct the transfer matrix, 
i.e., to consider a chain of the i?-operators and to take trace in the auxiliary space. We 
plan to address this problem elsewhere. It would be also very desirable to find a direct 
connection of our approach with elliptic beta integrals |2H[3l]. Presumably, the pairing 
(13. 2 p or (13. 4 p should be replaced by a sum of residues. 

Among other things, the results presented in this paper indicate convincingly that 
there should exist a meaningful theory of infinite-dimensional representations of the 
Sklyanin algebra. Such a theory is still to be developed and this paper may provide 
some background in reaching this ambitious goal. 

At last, one should keep in mind that the Sklyanin algebra is just a very particular 
representative of a wide family of elliptic algebras [36] and, moreover, integrable sys- 
tems associated to algebras from this class can be constructed |37]. It would be very 
interesting to investigate to what extent the methods developped in the present paper 
can be extended to other elliptic algebras and corresponding integrable models. Such an 
extension will probably require a further generalization of elliptic hypergeometric series. 
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Appendix A 

Theta-functions 

We use the following definition of the Jacobi ^-functions: 



feez 



e2iz\r) = E exp (-Kirik + + 2TTiz{k + - 
feez ^ ^ ^ 

93{z\t) = E exp (^irirk^ + 2Txizk 



64{z\r) = E exp ('nirk'^ + 2'ni[z H — )k 
fcez ^ ^ 



(Al) 



They also can be represented as infinite products. The infinite product representation 
for the 6*1(2;^) reads: 



00 

e,{z\r)=ie^-p[^-i'Kz) g(i_e2-'=^)(l-e2'^*«'=-i)^+^))(l-e2-(^^-^)). (A2) 



k=l 



Throughout the paper we write 0a(x|r) = 9a{x)^ ^(-^If) = ^{z)- The transformation 
properties for shifts by the periods are: 

ea{x ± 1) = (-l)^'"-^+^-2^„(x) , ea{x ± r) = (-l)^'"-i+^-4e-^^^^^^^^^„(a;) . (A3) 

Under the modular transformation r — t- —1/r the ^^-functions behave as follows: 

ei{z\T) = i y^e-^*^'/^^i(2/r| - 1/r) , 
e2{z\T) = yz7^e-^^^'/^^4(^/r| - 1/r) , 

(A4) 

93{z\t) = Jt/^e-^^^'/^9siz/r\ - 1/r) , 



e^izlr) = Ji/^e-^"^l^e2{z/T\ - 1/r) . 
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The identities often used in the computations are 

Hx%{y) + 0^{y%{x) = 2e,{x + y)9^{x - y), 

Ux)Uy) - OMUx) = 29,{x + y)9,{x - y), 

(A5) 

03ix%{y) + hivMx) = 293{x + y)9;{x - y), 

Ux)Uy) - Uy)U^) = 29^{x + y)9^{x - y), 
9i{z-a- d)9i{z -b- c)9i{a - d)9i{c - b) 
+ 9i{z - b - d))9i{z - a - c)9i{b - d)9i{a - c) (A6) 

= 9i{z-c- d)9i{z -a- b)9i{a - b)9i{c - d). 

By 6„ we denote the space of 6'-functions of order n, i.e., entire functions -F(x), x G C, 
such that 

F{x + l) = F{x) , F(x + r) = (-l)"e-™"-'™"F(a;) . (A7) 

It is easy to see that dim0„ = n. Let F{x) G On, then F{x) has a multiphcative 
representation of the form F{x) = cHILi^il^ ~ Xi), Yli=iXi = 0, where c is a constant. 
Imposing, in addition to flA7p . the condition F{—x) = F{x), we define the space 6^ C 9^ 
of even ^-functions of order n, which plays the important role in representations of the 
Sklyanin algebra. If n is an even number, then dimB^ = |n + 1. 



Elliptic gamma-function 

Here we collect the main formulas on the elliptic gamma-function [261 EZ] • We use the 
(slightly modified) notation of [27] . The elliptic gamma-function is defined by the double- 
infinite product 

oo ^ _ ^2ni{{k+l)T+(k'+l)T'-z) 

T{z\t,t')= n 1 _ .2.^(kr+k'r'+^) " (^^) 

k,k'=0 ^ ^ 

A sufficient condition for the product to be convergent is Imr > 0, Imr' > 0. We need 
the following properties of the elliptic gamma-function: 

T{z + l\T,r') = T{z\T,T'), (A9) 

T{z + t\t, t') = -ie~'^ri^\T')e''''9i{z\T')T{z\T, r') , (AlO) 
T{z + r'|r, r') = -ie~^r]o\T)e''''9i{z\T)T{z\T, r') , (All) 

where 

oo 

voir) = nil - e'-^'-) 

k=l 

is the Dedekind function. Another useful property is 

T(z\T,T')T{r'-z\r,r') = ——^^^. (A12) 
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Note also that T{z\r, t')T{t + r' - z\r, r') = 1. 

Under the modular transformation r — 1/t the elliptic gamma-function behaves 
as follows [27] : 



T{z/t\-1/t,t'/t) 

where 



T{z\t, r') = e-^(^) '^%V , , , , (A13) 

^ r((2-r)/r'| -r/r',-l/r') ^ ' 



1 3 r + r' - 1 2 r2 + r'2 + 3rr' - 3r - 3r' + 1 
^ ' 3rr' 2rr' 6rr' 

(AM) 

(r + r' - l)(r + r' - rr') 
12rr' ■ 

Let us list the most frequently used formulas for V{z) = r{z\T,2ri). Using (lAlip 
several times, we obtain: 

""^"".l ^^"^^ = e-^vk^R-k^-^k. TT ^^(^ ^ 2jr/) , (A15) 



r(x 2fcr/) ^ ._-^.k^^^r,k^J^k^-n^kx TT / g ^ 2r/ + 2jr])) \ (A16) 

where i? = In particular, ratios of such functions are expressed through 

the elliptic Pochhammer symbols as 

r{2ar] + 2kr]) _ n^i(^^B)kv ^i^^v) 



e 



r{2Pr] + 2kr]) r(2/3r/) 

(A17) 

r(2/3r7 - 2A;r7) r(2/3r7) [1 - a]^ ' 

As is seen from (lASp . the function r(2;|r, 2?]) has zeros at the points z = 2{k + I)// + 
(m + l)r + n, and simple poles at the points z = —2kri — mr + n, where k, m run over 
non-negative integers and n over all integers. The residues of the elliptic gamma-function 
at the poles at z = —2kri, k = 0,1,2, .. . are: 



k 

res _ r(^) = i-l)''e-'^^"R%f[(ei{2j7^)y' , (A18) 



z=—2kri 

where 

7ri(r+2r?)/12 



To = res 



^=0^^^^ 2mr]D{T)7]D{2r]) 



Elliptic hypergeometric series 

Here we follow [18j. We define the elliptic Pochhammer symbol (the shifted elliptic 
factorial) by 

[x],, = [x][x + l]...[x + k-l], (A19) 
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where [x] = 9i{2x7]) (cf. fl2.5p ). By definition, the elliptic hypergeometric series is 



This is an elliptic analog of the very- well-poised basic hypergeometric series [35]. The 
series is said to be balanced if z = 1 and 

r-2 

r - 5 + (r - 3)ai = 2 ^ a„+3 • (A21) 

m=l 



For a series I]a;>o c/c of the form f lA20p . the balancing condition f lA2ip means that the ratio 



Cfc+i/cfc of the coefficients is an elliptic function of k. For balanced series f lA20p . we drop 



the argument z = 1 and the parameters 77, r writing it simply as ,.+iCi;r-(ai; ^4, • • • , ctr+i)- 
For instance, 

The series is called terminating if at least one of the parameters a^, . . . , a^+i is equal 
to a negative integer number. In this case the sum is finite and there is no problem of 
convergence. If, say ctr+i = —n, then the series terminates at k = n. The terminating 
balanced series were shown [18] to possess nice modular properties. That is why they 
were called modular hypergeometric series. 

The modular hypergeometric series obey a number of impressive identities. One of 
them is the elliptic analog of the Jackson summation formula: 

[ai + l]„[ai-a4-a5 + l]n[ai-a4-a6 + l]n[ai-a5-a6 + l]n 
8Co'7(ai; ^4, . . . , 07, —n) = — — ^-^ — — — 

[ai -0:4 + Ij„[ai-a5 + lj„[ai-a6 + ljn[ai -04-05 -ae + ljn 

(A23) 

which is valid provided that the balancing condition 2ai + 1 = + + + — n is 
satisfied (the Frenkel-Turaev summation formula [IB]). 

A remark on the notation is in order. In the modern notation [19], what we call 
r+iujr{ai;a4, . . . ,ar+i\r],T) (following [18]), would be r+3K+i(ai; oe, • • • , flr+slo'^jP) with 
q = e^'^*'', p = e^'^*'^, aj = e^'^^'^"-^-'^ . In particular, our would be ^V^. We understand 
that the modern notation is better justified by the meaning of the elliptic very-well- 
poisedness condition than the old one and is really convenient in many cases. However, we 
decided to use the old Frenkel-Turaev notation for the reason that the additive parameters 
Q!j are more convenient for us than their exponentiated counterparts. We think that it 
is simpler than to introduce a version of r+iK with additive parameters. 



Appendix B 

In this appendix we give some details of calculations which involve modular hypergeo- 
metric series. 
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The normalization of the PV^-kernel 



Let us consider convolution of the kernels W^{X) and W^,{—\) given by equation (15.41) : 



- c(\)c( \) V f V ~ + + ^^^^ ^ 5(z z' + 2nn) 

In order to calculate it explicitly, consider the sum 

S(,)-T e,{2z-2X + 4kv) 

where the ly-functions are given by (I4.14p : 

y^,z-x+2kr,,zf. , -i^ix+v)(z-x+2kv) r(2z + 27] + 2kr,)V {2r^ + 2kr^) 

^ " r{2z-2X + 2k7])r{-2X + 2k7]y 



yyz+2nv,z-X+2kr,^^ _ A) = 6^ 



{X-ri){z+2nri) 



\'{2z-2X + 2r]+2n7]+2kr])r{2T] + 2nr] - 2k7]) 



Vi2z + 2n7] + 2kT])r{2X + 2n7] - 2kr]) 

Plugging this into ( IBip and representing ratios of elliptic gamma-functions through el- 
liptic Pochhammer symbols with the help of ( lAlSp . ( ]A16p . we obtain: 



SJz) = e 



4TTiz-^X{X+ri)-4TTi(X~ri)n 



^ g r(-2A)r(2z - 2A)r(2z + 2nr^)r(2A + 2nr^) 

^ ^ ^ r{2ri)r{2r] + 2nri)V(2z + 2ri)V{2z -2X + 2nr] + 2r]) 

^[^ + 2k][^], [-^]. [^+n]4-n]. 



X 



£ itMk [f + [^ + n + l],[-^^-n + l]. 



The sum in the last line is the terminating balanced elliptic hypergeometric series 

( z — X X z z — X + riz — X + ri \ 

suj ; , - + n, , , -n 

\ 7] 7] rj 2rj 2rj J 

which is equal to 

it + l]n [1 - n]„ [^]. [-^ - n + 1]„ 



[f + -n + lU^U-^ -n + ll 
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(see flA23l) ). Because of the factor [1 — n]„ this is zero unless n = 0. Therefore, Sn{z) = 
if > 1 and 



= ^^ni.-^-fx(x+v) r(2A)r(-2A)r(22:)r(22: - 2X)ei{2z - 2A) 



n{2r])r{2z + 2T])r{2z-2X + 27]) 
We thus have 

''dCW^{X)W^,{-\) = c{\)c{-\)9i{2z)So{z)S{z - z'). 

Using identities for the elhptic gamma-function the product 9i{2z)Sq{z) can be simphfied 
to 

e,{2z)S,{z) = Po ie-2-^^'/T(2A)r(-2A), 

where 

Po = — = TTT^- l^^J 



So, setting 



e(A) - ^ (B3) 



we obtain the relation /rfC (A)Vr5(-A) = 5{z - z'). 

The star-triangle relations 

Let us verify the star-triangle relation (17. 7p 

w'''^''{^l - v)w''^'"{\ - ^l)w':n{\ -v) = j dcw^ix - fi)w'''^^{\ - i^)w^„{ii - v) (B4) 

(see Fig. [TD]). We use formulas dHHD, ([53D- The left hand side is 
= c(A -iy)J2 Cn{z', z)6{z -z" -X + iy + 2nr]), 

n>0 

where c(A) is given by (1B3P and 

C„(z',z) = e-'fl(^-'^)^'-(^-^+^)^+(^-^)(^-'^+")l^i(2;2-2A + 2z/ + 4nr7) 

r(2i/ - 2A)r(2;z - 2A + 2iy)r{z+z'+fi-u+ri)r{z'-z+2X-i2-u+r]) 
r{2r]) V{2z + 2vi)V{z + z' - 2\ + ^i + v + r]) r{z' -z-fi + u + r]) 

[ ^ \n [ \n I 2ri J" I. 2r; J" 

Til [£ _l_ 11 2A+;^+/j+-f; i t z-z' -2\+v+^i+ri i 

[J-Jn -r IJra [ 2r? J" I. 2»7 J" 
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One can see from this expression that the left hand side of (IB4p is the kernel of the 
difference operator 

^2si(^x-u)iz-z')-^{\-uf \'i2z-2X + 2iy+2r])r{z+z'+fi-i'+r])r{z'-z+2X~fi-iy+ri) 

r{2z + 27]) r{z + z' - 2X + fi + u + T]) r{z' -z-^i + u + r]) 



e ^ 



V . ( z-^ + J^ . z + z' + u-fi + r] z-z' + u-fi + r] _ 



(B5) 



Let us turn to the right hand side of (IB4I) . It is 

w''^^{x-u)ei{2C)ei{2z") 



c(A — ij)c{iJ, — v) d( 



W<''{X - /i + 7])W'"''^{fI -U + 7]) 

X J2 S{z-C-X+fx + 2kri)S{C-z"-f^+iy+2k'ri) 

k,k'>0 

:{X - ]i)c{]i -v)Yl Bn{z', z) 5{z-z" -X + v+2n7]), 

n>0 



where 



" ei(2z-2X + 2iJi+ Akri)ei(2z-2X + 2v+ An'q)W^' ^^-^+^'+'^^'^{X-v) 

The next step is to identify this sum with the terminating elliptic hypergeometric se- 
ries with a pre-factor. The latter is essentially a product of ratios of the f-functions. 
Specifically, we have: 

^ g2s[(A-^)(2'-2)+(A-M)(A-M+»?)+(A-i'){A«-i'+r))]+47ri2:+47ri(M-i^+»?)n 



V{2ii-2X)V{2z-2X + 2ii)V{z + z'+ii-u+i])V{z'-z+2X-ii-u+i]) 
r(2?7) V{2z + 2ini)V{z+z'-2X-^ii+v+r])V{z'-z-ix + v+vi) 

. r(2^-2A+2.+2,,,) r(2.-2,+2n„) _ 
r(277 + 2nr/)r(2z-2A + 2^+2r/+2nr7) ^ 

X 8W7(ai;a4,---,a7, -'^), 

where the parameters are: 

Z — A + yU yU — A z — A + z/ 

«! = , a4 = , as = \- 

rj rj rj 

Z + z' + fl — U + T] Z — z' + fi — U + T] 

27] 27] 

The series with these parameters is balanced, so one can apply the Frenkel-Turaev sum- 
mation formula f lA23p . The result is 
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Now it is straightforward to calculate the ratio Cn{z' , z) / Bn{z' , z). One can see that all 
z, z' and n dependent factors cancel in the ratio and one is left with 

Cn{z', z) _ V{2r]) e^^^~^^^"~^^r(2;/-2A) _ c(A - /i)c(/i - z/) 
Bn{z',z)~ te^T^D{r) r(2/i - 2A)r(2z/ - 2/i) ~ c{\ - u) 

where c(A) is given by ( 1B3I) . This means that the left and right hand sides of (1B4I) are 
indeed equal to each other. 

The other star-triangle relation, (17. 8p . is proved in a similar way. We note that its 
both sides are kernels of the difference operator 



e "J 



V{2z-2\ + 2v+2'q)V{z+z'+\-ix+ri)V{z-z'+\-fi+if]) 



V{2z + 2?7) V{z + z' + 2v-\-^+ri) V{z-z'+2u-\-^+r]) 

[ z-\ + 



z-X + v _ u-X z + z'-X + fx + T] z-z'-X + fx + ri _ ^2r,a, j .^{u-m 
T] ^ rj ^ 2r] ' 2r] ) ' 



(B6) 



References 

[1] R.Baxter, Partition function of the eight-vertex lattice model, Ann. Phys. 70 (1972) 
193-228; 

R.Baxter, Eight-vertex model in lattice statistics and one- dimensional anisotropic 
Hezsenberg chain. I, II, III, Ann. Phys. 76 (1973) 1-24, 25-47, 48-71 

[2] L.D.Faddeev and L.A.Takhtadzhan, The quantum method of the inverse problem 
and the Heisenberg XYZ model (English translation), Russian Math. Surveys 34:5 
(1979) 11-68 

[3] T.Takebe, Generalized Bethe ansatz with the general spin representation of the 
Sklyamn algebra, J. Phys. A 25 (1992) 1071-1083; 

T.Takebe, Bethe ansatz for higher spin eight vertex models, J. Phys. A 28 (1995) 
6675-6706 

[4] E.K.Sklyanin, On some algebraic structures related to the Yang-Baxter equation. 
Funk. Anal, i ego Pril. 16:4 (1982) 27-34 

[5] E.K.Sklyanin, On some algebraic structures related to the Yang-Baxter equation. 
Representations of the quantum algebra. Funk. Anal i ego Pril. 17:4 (1983) 34-48 

[6] L.N.Lipatov, High energy asymptotics of multi-colour QCD and exactly solvable lat- 
tice models, JETP Lett. 59 (1994) 571-574 

[7] L.Faddeev and G.Korchemsky, High energy QCD as a completely integrable model, 
Phys. Lett. B342 (1995) 311-322 

[8] S.Derkachov, G.Korchemsky and A.Manashov, Noncompact Heisenberg spin mag- 
nets from high energy QCD, I: Baxter Q- operator and separation of variables, Nucl. 



30 



Phys. B617 (2001) 375-440; 

M.Kirch and A.Manashov, Noncompact SL{2, R) spin chain, JHEP 0406 (2004) 
035; 

S.Derkachov, Factorization of the R-matrix and Baxter's Q-operator, J. Math. Sci. 
151 (2008) 2848-2858, arXiv:math.QA/0507252, 

[9] S.Derkachov, D.Kharakhanyan and R.Kirschner, Baxter Q-operator of the XXZ 
chain and R-matrix factorization, Nucl. Phys. B738 (2006) 368-390; 
A.Bytsko and J.Teschner, Quantization of models with non-compact quantum group 
symmetry. Modular XXZ magnet and lattice sinh- Gordon model, J. Phys. A: Math. 
Gen. 39 (2006) 12927-12981 

[10] I.G.Korepanov, Vacuum curves of C-operators connected with the six-vertex model 
and construction of TZ-operators, Chelyabinsk Polytechn. Inst., Deposited in the 
VINITI 02.04.86, N 2271-B86; Algebra i Analiz 6:2 (1994) 176-194 

[11] V.V.Bazhanov and Yu.G.Stroganov, Chiral Potts model as a descendant of the six- 
vertex model, J. Stat. Phys. 59 (1990) 799-817 

[12] R.Kashaev, V.Mangazeev and T.Nakanishi, Yang-Baxter equation for the sl{n) chi- 
ral Potts model, Nucl. Phys. B362 (1991) 563-582 

[13] K.Hasegawa and Y.Yamada, Algebraic derivation of the broken Zj^i- symmetric model, 
Phys. Lett. A 146 (1990) 387-396 

[14] Y.-H.Quano and A.Fujii, Yang-Baxter equation for broken Z^*^"^ models. Mod. Phys. 
Lett. A 17 (1993) 1585-1597; 

Y.-H.Quano, Generalized Sklyanin algebra and integrable latice models. Int. J. Mod. 
Phys. A 9 (1994) 2245-2281 

[15] K.Hasegawa, L-operator for Belavin's R-matrix acting on the space of theta func- 
tions, J. Math. Phys. 35 (1994) 6158-6171 

[16] A.Kirillov and N.Reshetikhin, Representations of the algebra Uq{sl2), q-orthogonal 
polynomials and invariants of links, in: Infinite dimensional Lie algebras and groups, 
ed. V.Kac, 285-339, Adv. Ser. in Math. Phys. 7, World Scientific, Singapore, 1988 

[17] A.Zabrodin, Commuting difference operators with elliptic coefficients from Baxter's 
vacuum vectors, J. Phys. A: Math. Gen. 33 (2000) 3825-3850 

[18] I.Frenkel and V.Turaev, Elliptic solutions of the Yang-Baxter equation and mod- 
ular hypergeometric functions, in: The Arnold-Gelfand Mathematical Seminars, 
Birkhauser, Boston, MA, 1997, 171-204 

[19] V.Spiridonov, Classical elliptic hypergeometric functions and their applications, 
Rokko Lect. in Math., Vol. 18, Dept. of Math, Kobe Univ., 2005, pp. 253-287; 
V.Spiridonov, Essays on the theory of elliptic hypergeometric functions, Uspekhi 
Mat. Nauk 63 no 3 (2008) 3-72 (Russian Math. Surveys 63 no. 3 (2008) 405-472); 
V.Spiridonov, Continuous biorthogonality of the elliptic hypergeometric function. Al- 
gebra i Analiz (St. Petersburg Math. J.) 20 (2008) 155-185 



31 



[20] S. Ole Warnaar, Summation and transformation formulas for elliptic hypergeometric 
series, Constructive Approximation 18 (2002) 479-502; 

S. Ole Warnaar, Summation formulae for elliptic hypergeometric series, Proc. Amer. 
Math. Soc. 133 (2005) 519-527 

H. Rosengren, An elementary approach to 6j-symbols (classical, quantum, rational, 
trigonometric, and elliptic), Ramanujan J. 13 (2007) 133-168 

S.Derkachov, D.Karakhanyan and R.Kirschner, Yang-Baxter R operators and pa- 
rameter permutations, Nucl. Phys. B785 (2007) 263-285, |arXrv:hep-th/0703076| 

V.Bazhanov and S.Sergeev, A master solution of the quantum Yang-Baxter equation 
and classical discrete integrable equations, larXiv: 1006.06511 

V.Spiridonov, Elliptic beta integrals and solvable models of statistical mechanics, 
'arXiv:1011.3798 

A.Zabrodin, On the spectral curve of the difference Lame operator. Int. Math. Res, 
Not. 1999 (1999) 589-614 

S.N.M.Ruijsenaars, First order analytic difference equations and integrable quantum 
systems, J. Math. Phys. 38 (1997) 1069-1146 

G. Felder and A.Varchenko, The elliptic gamma function and SL{3, Z) x Z^, Adv. 
Math. 156 (2000) 44-76 

I. Krichever, Baxter's equations and algebraic geometry. Funk. Anal, i ego Pril. 15:2 
(1981) 22-35 

I.Krichever and A.Zabrodin, Spin generalization of the Ruijsenaars-Schneider model, 
non-abelian 2D Toda chain and representations of Sklyanin algebra, Usp. Mat. Nauk, 
50:6 (1995) 3-56 

I.Krichever and A.Zabrodin, Vacuum curves of elliptic L-operators and representa- 
tions of the Sklyanin algebra, Amer. Math. Soc. Transl. (2) 191 (1999) 199-221 

H. Rosengren, Sklyanin invariant integration. Int. Math. Res. Not. 2004 (2004), 3207- 
3232 

E.Rains, Transformations of elliptic hypergeometric integrals, Ann. of Math. 171 
(2010) 169-243, arXiv:math.QA/0309252 

E.Date, M.Jimbo, T.Miwa and M.Okado, Fusion of the eight-vertex SOS model, 
Lett. Math. Phys. 12 (1986) 209-215, Erratum and Addendum: Lett. Math. Phys. 
14 (1987) 97; 

E.Date, M.Jimbo, A.Kuniba, T.Miwa and M.Okado, Exactly solvable SOS models, 
II: Proof of the star-triangle relation and combinatorial identities. Advanced Studies 
in Pure Math. 16 (1988) 17-122 

[34] V.Spiridonov, On the elliptic beta function, Uspekhi Mat. Nauk 56 (2001) 181-182 
(Russian Math. Surveys 56 (2001) 185-186) 



32 



[35] G. Gasper and M.Rahman, Basic hypergeometric series, Encyclopedia of Mathemat- 
ics and its Apphcations, Volume 35, Cambridge University Press, 1990 

[36] A.Odesskii and B.Feigin, Sklyanin elliptic algebras, Funct. Anal. Appl. 23 (1989) 
207-214; 

A.Odesskii, Elliptic algebras, Russian Math. Surveys 57 (2002) 1127-1162 

[37] A.Odesskii and V.Rubtsov, Integrable systems associated with elliptic algebras, 
Rokko Lectures in Mathematics 18 (2005) 111-132, eds. M.Noumi and K.Takasaki, 
[arXiv:math /0404159 



33 



